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Abstract 


In  this  report  the  aerodynamic  pressure  acting  on  a  circular 
cylindrical  shell  of  infinite  length  in  an  airflow  parallel  to  the  cylin¬ 
der  axis  is  studied.  The  shell  is  deformed  by  a  harmonically  oscil¬ 
lating  standing  wave  of  a  sinusoidal  pattern.  Based  upon  the  exact 
solutions ,  asymptotic  expansions  are  developed  for  the  aerodynamic 
pressure.  In  this  manner  the  accuracy  of  the  linear  piston  theory 
approximation,  when  applied  to  cylindrical  shells,  is  investigated. 
Furthermore  improved  approximations  can  be  obtained  from  these 
asymptotic  expansions. 
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Velocity  of  sound  in  the  undisturbed  air. 
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Radial  coordinate. 

sin(vx)  or  cos(vx),  see  eq.  (2.4b). 

Time. 

Deformation  of  the  shell,  positive  in  the 
positive  r -direction. 

Real  or  complex  constant. 

See  eq.  (6.2). 

Length  coordinate. 

Cylinder  coordinate  system. 

Angular  coordinate. 

Wave  number,  see  eq.  (1.4). 
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Density  of  the  undisturbed  air. 

Velocity  potential. 

Circular  frequency  of  the  shell  vibration, 
see  eq.  (2.1),  real  and  non-negative. 
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INTRODUCTION 

Piston  theory,  originally  developed  by  Lighthill  (ref.  10), 
was  introduced  into  aeroelasticity  in  the  linearized  form  by  Ashley 
and  Zartarian  as  a  handy  tool  in  1956,  see  ref.  2.  This  theory 
furnishes  an  approximation  for  the  aerodynamic  pressure  acting 
on  a  slightly  deformed  flat  plate  in  a  supersonic  air  stream.  The 
linearised  piston  theory  is  widely  used  in  the  investigation  of  the 
flutter  of  flat  panels,  see  the  summaries  of  Fung,  refs.  5  and  6. 

As  far  as  the  flutter  of  cylindrical  shells  is  concerned,  an 
appropriate  approximation  for  the  aerodynamic  pressure  acting 
on  a  vibrating  shell  is  lacking.  In  the  case  of  supersonic  outer 
air  stream,  many  authors  use  the  linear  piston  theory  expression 
for  the  aerodynamic  pressure.  Some  authors  even  employ  Ackeret's 
formula,  neglecting  the  aerodynamic  damping  term  of  the  linear 
piston  theory. 

There  are  doubts  about  the  accuracy  of  using  the  linear 
piston  theory  for  cylindrical  shells.  This  opinion  was  strengthened 
by  recent  studies  at  Caltech.  The  author  has  investigated  the  flut¬ 
ter  in  an  axisymmetrical  mode  of  a  thin  cylindrical  shell  of  finite 
length  in  a  supersonic  airstream  on  the  bases  of  Timoshenko's 
linearized  shell  equations,  linear  piston  theory  and  linear  material 
damping,  see  ref.  7.  On  the  other  hand,  wind  tunnel  experiments 
in  the  summers  of  1961  and  1962  (see  ref.  11)  reveal  no  flutter  in  an 
axisymmetrical  mode  in  the  predicted  Mach  number  range.  This 
indicates  that  the  basic  physical  assumptions  on  which  the  above 
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mentioned  calculations  are  based  are  not  satisfactory.  The  inade¬ 
quacy  of  the  linear  piston  theory  and  the  neglect  of  boundary  layer 
effects  might  be  responsible  at  least  partially  for  this  disagreement. 

The  influence  of  the  boundary  layer,  using  a  simplified  model 
set  forth  by  Fung  in  the  summer  of  1961,  has  been  studied  by  Fung 
and  Anderson,  see  refs.  1  and  6.  It  is  the  subject  of  the  present 
paper  to  investigate  the  accuracy  of  the  linear  piston  theory  approx¬ 
imation  when  applied  to  cylindrical  shells. 

Leonhard  and  Hedgepeth  develop  in  ref.  9  an  exact  expres¬ 
sion  for  the  aerodynamic  pressure,  acting  on  a  cylindrical  shell  of 
infinite  length,  which  is  exposed  externally  to  an  air  stream  parallel 
to  the  generators  of  the  cylinder,  and  where  the  shell  is  slightly 
deformed  by  a  harmonically  oscillating  standing  sinusoidal  wave. 

In  the  present  paper  the  well-known  asymptotic  expansions  for  cyl¬ 
inder  functions  are  used  to  obtain  an  asymptotic  formula  for  the 
aerodynamic  pressure  under  consideration,  in  the  form  of  poly¬ 
nomials  of  the  reciprocal  of  the  shell  radius  with  a  remainder  term. 
The  results  differ  according  to  whether  the  Mach  number  of  the 
free  stream  or  the  frequency  of  the  shell  vibration  is  "large"  or  not. 

In  the  case  of  "large"  Mach  numbers  or  "large"  frequencies, 
the  asymptotic  expansion  can  be  Bplit  up  into  the  linear  piston 
theory  expression  and  some  correction  terms.  The  latter  ones 
tend  to  zero  if  the  radius  of  the  shell  and  the  Mach  number  tend  to 
infinity;  the  same  is  true  if  the  shell  radius  and  the  frequency  of  the 
shell  vibration  tend  to  infinity.  Hence  this  expansion  provides  an 
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estimate  of  the  error  of  the  linear  piston  theory  approximation. 
Furthermore  improved  approximations  can  be  obtained  from  this 
asymptotic  expansion. 

In  case  neither  the  Mach  number  nor  the  frequency  is  "large" , 
the  linear  piston  theory  approximation  can  no  longer  be  considered 
as  a  first  order  approximation  of  the  considered  aerodynamic  pres¬ 
sure. 

Furthermore,  on  the  assumption  that  the  Mach  number  is 
larger  than  one  and  letting  the  shell  radius  tend  to  infinity  and  the 
frequency  of  the  standing  wave  tend  to  zero,  the  asymptotic  expan¬ 
sion  tends  to  Ackeret's  formula. 

These  investigations  have  been  performed  for  standing 
waves  of  a  sinusoidal  pattern  on  the  cylinder.  The  exact  terms 
for  the  corresponding  aerodynamic  pressure,  as  well  as  most  of 
the  coefficients  of  the  asymptotic  expansion,  depend  also  upon  the 
frequency  of  the  shell  vibration  and  upon  the  wave  length  in  the 
direction  of  the  cylinder  axis.  Therefore,  if  these  results  are  to 
be  applied  to  more  general  vibration  patterns,  every  term  of  the 
corresponding  Fourier -expans ion  has  to  be  treated  specially. 

It  is  with  pleasure  that  the  author  expresses  his  gratitude 
to  Professor  Y.  C.  Fung  for  fruitful  discussions  and  for  bringing 
the  above-mentioned  paper  of  Leonhard  and  Hedgepeth  again  to  the 
attention  of  the  author.  Furthermore  the  author  wishes  to  express 
his  thanks  to  Mrs.  O.  M.  Eaton  for  her  careful  numerical  calcu¬ 
lations  and  to  Mrs.  E.  Fox,  who  again  typed  the  manuscript. 
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1.  Aerodynamic  Pressure  Acting  on  a  Stationary  Sinusoidally  De¬ 
formed  Cylindrical  Shell 

We  consider  an  infinitely  long  circular  cylindrical  shell  of 
radius  R,  which  is  exposed  externally  to  a  uniform  air  stream 
parallel  to  the  generators  of  the  cylinder.  The  Mach  number  of  the 
undisturbed  airstream  is  M.  The  density  and  the  velocity  of  sound 
of  the  air  are  denoted  by  po  and  ao  respectively.  Let  a  cylinder 
coordinate  system  x,  r,  8  be  chosen,  where  the  positive  direction 
of  the  x-axis  coincides  with  the  positive  direction  of  the  airstream^. 
Now  the  shell  is  assumed  to  be  stationary  and  slightly  deformed  by 
a  sinusoidal  wave 

w  ■  w0  cos  (n0)  e"*vx  (1.1) 

where  w  is  the  radial  displacement  of  the  shell,  measured  positive 
in  the  outside  direction.  w0  is  a  real  or  complex  number,  n  is  a 
non-negative  integer  and  v  is  a  positive  or  negative  number.  The 
deformation  (1.1)  of  the  shell  is  assumed  to  be  of  infinitesimal  ampli¬ 
tude,  so  that  the  differential  equations  and  the  boundary  conditions 
for  the  air -flow  are  linearized,  and  the  boundary  conditions  are  ap¬ 
plied  on  the  mean  position  of  the  shell  rather  than  on  the  actual  sur¬ 
face.  The  governing  equation  for  the  velocity  potential  4>{x,r,0) 
reads  (see  e.g.  ref.  4,  page  432)  for  Jm  j  #  1: 

^  Therefore  a  negative  Mach  number  means  that  the  airstream  is 
moving  in  the  negative  x-direction. 
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(i-M*, »!$ t  *  i  * 0.2) 

Bx2  dr2  r  9r  x2  802 
The  boundary  conditions  are  given  by 

for  r  ■  R; 

(1.3) 

and  Sommerfeld' s  finiteness  and  radiation  condition  for  r-*co. 


Solutions  of  eq.  (1.  2)  which  are  suitable  for  our  purposes  are  of  the 
form 


-1VX 


0(x,r,8)  *  f (r )  cos  (n0)  e 


where  f(r)  is  a  solution  of  the  differential  equation 

^  +  7  3!  +  (vW-l)  -  f(*>  =  0 

dr  r  dr  r 


a.  4) 


(1.5) 


and  satisfies  boundary  conditions  corresponding  to  the  conditions  (1.3). 
Knowing  the  velocity  potential  (f>  one  obtains  the  resulting  pressure 
perturbation  Ap(x,r,8)  by  (see  e.g.  ref.  4,  page  432) 


Ap(x,r,6)  «  -  p0  a,,  M  . 


a.6) 


In  this  manner  one  arrives  at  the  following  results  for  the  pressure 
perturbation  acting  on  the  cylinder  well,  Ap  ■  Ap(x,R,0),  which  de¬ 
pends  upon  the  Mach  number  range  under  consideration  (for  details 


the  reader  is  referred  to  ref .  9): 

Po(a0M)2  2  K  (|  vjR|/l-M2  ) 

*P  «  ;■  7-  °  L  _ —  w(x, 8) 

1 1-M2  Ivl  K^(  J  v  |R  fl-M2) 

for  -1  <  M  <  1  , 


(1) 


(1.7a) 


Try 


Here  and  in  the  following  we  always  choose  the  positive  square -root. 


w(x,e) 


-6- 

\l Mz-1  M  Jr  \^i2-i) 


for  M|  >  1  and  Mv  >  0, 


for  M  >  1  and  Mv  <  0  . 


(1.7b) 


ft.  7c) 


Here  H^(z)  and  H^\z)  are  the  Hankel  functions  (Bessel  functions 
of  the  third  kind)  and  Kn(z)  the  modified  Hankel  functions  (modified 
Bessel  functions  of  the  third  kind).  Because  Kn(z)  is  real-valued 
when  n  is  real  and  z  is  positive  (see  e.g.  ref.  3),  which  is  obvi¬ 
ously  the  case  here,  the  perturbation  pressure  Ap  for  |  M|  <  1  is 
“in  phase"  with  the  shell  deformation  w(x,0),  see  eq.  (1.7a).  This 
is  no  longer  true  for  J  M  J  >  1.  For  the  fractions  of  the  Hankel-functions 
and  their  derivatives  involved  in  eqs.  (1.7b)  and  (1.7c)  the  following 
relation  holds 


(1.8) 


for  j  M  |  >  1.  ^ 


^  Here  and  in  the  following  the  barred  symbol  denotes  the  conjugate 
complex  number,  i.e.  a+ib  *=  a-ib;  (a,b  real). 
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2.  Aerodynamic  Pressure  Acting  on  a  Shell  which  is  Deformed  by 
a  Traveling  Sinusoidal  Wave  or  by  a  Standing  Sinaeoidal  Wave 

Now  we  consider  shell  deformation  of  the  form 


w(x,6,t)  ■  w0  cos  (n6)  e’*vx  e*wt  * 


-iv(x  -£t) 


e  w0  cos  (n6)  e 


(2.1) 


.  («  >  0) 


where  u  is  assumed  to  be  a  real,  non-negative  number.  The  defor¬ 
mation  w  represents  a  sinusoidal  wave  which  travels  with  the 
velocity  ^  in  the  positive  x-direction.  A  moving  observer  who 
travels  with  the  same  velocity  "  in  the  positive  x-direction  can 
apply  the  formulas  (1.  7)  after  M  has  been  replaced  by 


Mi 


*  M  - 


OJ 

va0 


(2.2) 


Referring  these  results  again  to  the  (resting)  x,  r, 6-coordinate  sys¬ 
tem  one  obtains  for  the  pressure  perturbation  acting  on  the  vibrating 
cylinder  wall 


Ap 


Po  (ao  M^)2  y2 
^l-Mj2  M 


-iv(x-2t) 

w0  cos(n8)e 


for  |mJ<1, 


(2.3a) 


Ap 


PotaoMj)2  y2  Hjl)(|v  Ir/m2-!) 

)/m2  -  1  lv 


«(!)'( | 


w0  cos(n6)e 


-iv(x-St) 


for  J  Kdj  J  >  1  and  MjV  >  0  , 


(2.3b) 
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PolaoM/  v2 

h 


-iv(x-2t) 

w0  cos(n0)e 


(2.3c) 


for  |  Mj  >  1  and  MjV  <  0  . 

Shell  deformations  of  the  form 

w(x,0,t)  *  w0  cos  (n$)  »(*y  «iWt  <^-4*4 


with 


s(x)  ■ 


rsin  (vx) 
or 

cos  (vx)  , 


(2.4b) 


which  represent  a  standing  sinusoidal  wave,  are  obtained  by  suitable 
superpositions  of  two  waves  of  the  form  (2.1),  traveling  in  opposite 
directions.  Hence  the  corresponding  aerodynamic  pressure  is  a  su¬ 
perposition  of  the  corresponding  pressure  terms  (2.3).  Introducing 
the  notation 


M2  *  M  + 


6? 

vao 


(2.5) 


we  find  especially; 


Ap  *  iwo  cos(n 0)e 


icot 


_m|_  gwR* 


e 


-ivx 


(2.6a) 
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for 


Mx  >  1; 


M2  >  1;  v 


>  0;  w  *  w0  cos(n6)  sin(vx)e 


i<ot 


Ap  ■  iwo  co»(n8)e 


iwt 


„  2 
Poao 


(2.6b) 


(2.7a) 


for 


>  1;  M2 


>  1;  v  <  0;  w  *  w0  cos(n0)sin(vx)e 


icot^ 


(2.7b) 


(2.8a) 


for 


Mj  >  1;  M2  >  1;  v  >  0;  w  «  w0  coa(n0)cos(vx)e 


iut 


*  ,  iut  2  v 

Ap  ■  iw0  cos(n0)e  Po  a* 


;m l  h<2)(|v|r/^T)  ,iv. 


(2. 8b) 


(2.9a) 


-10- 


for 


iwt. 


<  -1;  M2  >  li  v  >  0;  w  *  w0  cos(n8)sin(vx)e  ; 

t4  e-iv* + 

2  V*i 1  Hi2>' '  I |/M? 'l) 


(2.9b) 


2  H  f  3 

Ap  •  wo  co«(n8)e  Po  *0  “7» —  i  i— — 


(2.10a) 


m2  42)Uv  1rVm|-d 

^ M2-l  HfdviR^l) 


for 


iwt. 


Mi 


>  1;  M  <  -1;  v  <  0;  w  *  wo  cos(n0)cos(vx)e  ; 


(2.10b) 


iwt  2  I  v  !  f  ^1 
Ap  «  iwo  co»(n8)e  Po^o  ^  ■<  — 


f  KJ|»  \Rjjl-M?) 


-ivx 


(2.11a) 


for 


j  j*!  |  <  1;  j  M2  |  <  l;  v  2  0;  w  ■  w0  cog(ne)  sin(vx)e  ; 

2 .  r jl.  VjlWES: 

Ap  *  iwo  cos(n8)«  Po&«  'J  |  F  I  V 

K;(|v|Rjl-M,2] 


iwt;  (2 .  lib) 


-IVX 


m2  h|i2)(  1  v  |r/m2-i)  #ivx| 
fuf 1  HfdvjR/^l, 


(2.12») 
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for 

J  Mjj  <  1;  >  1.  v  >  0;  w  *  w0  cos(n0)  sin(vx)  e1*1*  .  (2.12b) 

Letting  the  Mach  number  M  tend  to  zero,  one  obtains  the 
aerodynamic  pressure  acting  on  the  vibrating  cylinder,  which  per¬ 
forms  the  vibration  (2.4a)  in  still  air.  This  can  be  realised  by  in¬ 
spection  of  the  formulas  (2.9)  and  (2.11)  and  by  a  comparison  with 
the  corresponding  results  in  ref.  8.  (Notice:  The  positive  w-direc- 
tion  in  ref.  8  is  opposite  to  the  one  chosen  in  the  present  paper.) 
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3.  Asymptotic  Expansion!  of  the  Involved  Cylinder  Function  Term* 

In  this  section  we  apply  the  following  well-known  asymptotic 
expansions  for  cylinder  functions  for  4  —  oo  (see  e.g.  ref.  3,  pages 
85  and  86). 

{T  ~^r+0(e’Q)}.  (3.1) 

L  Vo  ('2iS)  J 


i(€  -  5jr-  3)  r?''1  (n*<l) 


Knte) 


” -4 
■?fe 


{Y  -^4+  o(rQ)}  . 


with 


(n.O)  =  1; 


(n.q)  «  (Z2\\)  (4n2-l)  (4n2-32)  ...  (4n2- [2q-l]  2); 


(3.2) 


(3.3) 


to  the  exact  expressions  for  the  aerodynamic  pressure  stated  in  sec¬ 
tion  2.  Using  eq.  (3.1)  one  arrives  in  a  straightforward  manner  at 


H^d) 

H^1}l(e) 


it 


n 


u-2 


/  2  3.  e-3  . 

(n  -¥)  4  + 


+0(|-4)  for  +oo. 


(3.4) 


Because  |  is  assumed  to  be  real  valued,  one  obtains  (see  eq.  (1.8)) 

2 


h£2)(4) 


(4) 


(  n  .  j  ,-l  ,  .  ,n  3,  t-2  ,_2  3t  ,-3  . 

vpr iS  nlT-5>«  + 


n 


(3.5) 


+  0(4"  )  for  |  ao 


^Here  and  in  the  following  it  is  assumed  that  %  is  a  real  variable.  O 
is  the  Landau -symbol;  i.e.  F(£)  ■  0(4~c*)  for  4  -*■«>  means  that 
F(4)4a  is  bounded  as  4  tends  to  infinity. 
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Further,  applying  eq.  (3.2)  one  finds 


*  -i  +  i  e_1  +  (^- 1)  e-2  -  (n2  - 1)  r3  +  ou*4) 
Kk(S> 

for  4  —  oo  . 


From  the  well-known  relations  (see  e.g.  ref.  3) 

H^g)  -  {-!)■»*>  (g);  H*2>(g)  (-l^^g); 

Kn(g)  -  K.n(g) 


(3.7) 


one  learns  the  following:  If  one  includes  in  the  asymptotic  expansions 
(3.4).  (3.5)  and  (3.6)  higher  powers  of  ,  then  the  corresponding 
coefficients  depend  on  powers  of  n  with  even  exponents. 

Table  I  on  page  28  gives  some  information  about  the  accuracy 
of  the  expansion  (3.4).  The  parameter  values  R  and  v  are  chosen  with 
respect  to  an  electroplated  copper  cylinder  which  was  used  repeatedly 
during  the  flutter  experiments,  see  ref.  11.  The  length  of  this  cylin¬ 
der  is  L  *  16  inches,  its  radius  is  Ri8  inches. 

Using  the  expansions  (3.4)  and  (3.5)  one  realizes  immediately 
that  for  R  -*  ao  the  aerodynamic  pressures  (1.7b)  and  (1.  7c)  (i.e. 
the  aerodynamic  pressure  acting  on  a  stationary  and  sinusoidally  de¬ 
formed  cylinder  for  |  M  |  >1)  tends  to  Ackeret*  s  formula 
*  Po(aoM)2  d(w0e'ivx) 

AP*  «  - - -  .  (3.8) 

/m2-i  <** 
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4.  Asymptotic  Expansion  o £  the  Aerodynamic  Pressure  lor  |  M^|>  1 


and  M,  ]  >  1. 


Let  the  cylinder  deformation  w  be  given  by  eq.  (2.4a)  and 
let  and  M2  be  restricted  by  J  |  >  1  and  |M2j>l«  We  apply 
the  results  of  section  3  to  the  exact  expressions  for  the  aerodynamic 
pressure,  which  are  partially  listed  in  eqs.  (2.6)  to  (2.10).  We  obtain 

Ap  «  ao  Po  {ao  M  + 

M,  M-, 


x  2  f.  /  1  ,  2  2<o\„.  L/  ,  2  ow\®wx 

‘-'•{’f-  + ‘  vtzit  * 

hr  m. 


— 

M2 


^ (■  rX" - T^r) w  *  SR 


— 7 

Mi 


— 7 
M2 


■  '■  l  —  '—n 


8w 


n2-3 

+  i_ 31 - Mn  1  ]V*. 


( _ _L_ 

2vR*  V  M,[l--iy]' 

“T 

2  3 
n  -? 


1 


I  w  + 


2l 

M2 


2  3 
n  -f 

~ 2T3 


M1 

M 

’  1 

1 

v4Rj  v  [1  *  "41  M2^  • 


m: 


m: 


2  3 


w  + 


aw, 


..  n  ‘I  /  1  .  1  \ 

m“[i- 


m; 


m: 


(4.1a) 


twK^  mJu  -  v  ^  )  *  o(^,m?u-  Xi572^' 

Ml 


i|[l-  -SI’ 


Mj  M2 


I 
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with 


Mz-M+7ir 


(4.1b) 


and 


w  m  w0  cos(n0)a(x)el“t;  s(x)  * 


sin(vx) 
or  0 

coa(vx) 


(4.1c) 


for 


Mj  |  >  1;  |  M2  J  >  1. 


(4.  Id) 


(Notice:  For  <  -1  it  holds  M^/[M^-1]^  “  -Mj/ll  -  — )  The 
first  term  on  the  right  hand  side  of  eq.  (4.1a)  is  the  approximation 
by  means  of  the  linear  piston  theory 


Ap  «  aop  o  -|ao  M  ^  .  (4. 2) 

which  has  been  separated.  The  other  terms  on  the  right  side  of  eq. 
(4.1a)  are  the  "correction  terms"  and  a  "remainder",  which  demon¬ 
strate  how  far  the  linear  piston  theory  approximation  (4. 2)  is  in  error. 
The  first  two  correction  terms  are 


1  2 

tao  Pqiv 


(- 


M. 


[1-^4-]* 


[i- 


m2 

~4rl 

MI 


i 


2  w  \ 
vao  / 


w«Fo  w 


(4.3a) 


i 
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The  remaining  correction  terms  are  of  the  form 


iao  Po 


m 

m -Inin 


“M1"' 

m2 


[1-  -^1  ’ 


(4. 3c) 


1  2 
2a0  PO 


vmRm  V^-lrj. 


M2 


■)£• 


(4.3d) 


m  =  1,  2,  3. 


m  9w 
,m  ~Zk 


*  2 
The  coefficients  and  Am  ,  m  =  2,  3,  are  polynomials  in  n  , 

see  section  3.  By  means  of  eqs.  (4.3)  new  coefficients  and 

jjt 

Fm  ,  m  *  0,  1,2,3,  are  introduced.  Naturally,  always  only  one  of 

the  two  signs  in  (4.3c)  and  (4.3d)  is  relevant.  Expanding  the  coeffi- 

*  11 
cients  Fm  and  Fm  in  power  series  of  and  -jj-  one  obtains, 

after  some  regrouping 


i  M|-Mf 


(4.4.) 


)  I  denotes  absolute  value  of  the  binomial-coefficient 
The  same  holds  for  the  following. 
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4b) 


Fj  * 


■5  ao  Po  H 


M2+M*  mJ+m* 


2  +  “TX  + 


2+  i  +  1 

^+***-j  ■ 


M’  M’ 

i  2.  f  «*7ST,  M2-Ml4  .  M2-“l6  .  1 

,2,  Am 


(4.4c) 


(4-4d) 


F  « 
m 


(4. 4e) 


,/- ^ <  3i< 3 .  •, 

*  2  )'  J 


m 


l2  A.  r^'-K"'1  . 

*  v“  l  ''  1  '* 


m+1  .  .m+  3 


033F-} ' 


(4.4£) 


Furthermore  it  holds 


(4.5a) 


Mjt  Mj  -  z{m*  +  Q  M**'2^)  +  .  .  .}  , 
“z-K*2^  {0  *“*  *  (3)  3 <5^  +  •  •  • } 

<“2  *  {m2  '  l-Tfe'2}1* ' 


(4.5c) 


ji  ■  1,2 
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where,  of  course,  the  dots  in  (4.5a)  and  (4.5b)  stand  for  finite  sums. 

Using  the  coefficients  F  and  F_  ,  introduced  by  eqs.  (4.3), 

m  m 

one  can  rewrite  eq.  (4.1a)  in  the  form 

Ap  b  aoPo /ao  M  -gjj  +  +  (4.6) 


+  F0w  +  F0*  +  £  (fx  w  +  F*  |~  )+ 


4-  w  1 

r\  (  1  1  \ 

-L/->  f  ^  1  1  I  X 

T  W  1 

-  '  v3R4  M*[l  -  -V-l5/Z' 

1  K 

'TiT*  m*[i  - 

M2 

dw 

rn  /  i  i 

L  'v4R4  M?[l  -  -L] 5/2 
M1 

M2 

Now  we  are  ready  to  discuss  the  behavior  of  the  correction 
terms  and  the  remainder  as  one  of  the  parameters  |mJ,  R,  <■>,  and 
jv  J  tend  to  -t-  oo,  provided  that  the  inequalities  (4.  Id)  are  observed. 
First  of  all  an  inspection  of  the  remainder  in  eq.  (4.1a)  or  eq.  (4.6) 
discloses  immediately  that  the  remainder  tends  always  to  zero  as 
soon  as  at  least  one  of  the  above  mentioned  four  parameters  tends 
to  +  oo. 

The  following  statements  are  based  on  eqs.  (4.1)  to  (4.6)  and 
the  just  quoted  statement  concerning  the  behavior  of  the  remainder. 
For  j  M  |  -»  oo  we  obtain 


F^-*  0  for  m  «  0, 2, 3;  F^  -*  0  for  m*  0,1, 2, 3; 


-  i»«Po  tor  j  M  j  —  oo, 


(4.7a) 
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hence 


Ap  -  {^!*  for  jM  j  -•>  oo. 


(4.7b) 


For  R  -*  oo  it  holds 
* 


F  F 

“  —  0;  —  —  0  for  R  —  oo,  m  ■  1,2,3, 
Rm 


TT 

m 


hence 


'0'^* 

•i6i^--}s- 


m^m2 


(4.8a) 


(4.8b) 


for  R  -*  oo. 


From  eqs.  (4.7)  and  (4.8)  we  conclude 

i  I 

Ap  -  Ap  -►  0  for  |M  j  -*  oo,  R  -*■  co. 


(4.9) 


Further  we  obtain 


F 


m 


0,  m  ■  0, 2, 3; 


0,  m  =  0,1, 2, 3;  Fj 


-  i«oPo 

(4.10a) 


for  «  —  +  oo , 


hence 
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2 

Ap  -  {Ap*  -  w}  -*  0  for  (i)  +oo.  (4.10b) 

Eqs.  (4.8)  and  (4.10)  lead  to 

$  l|( 

Ap  -  Ap  _  0  for  R  — ►  cd,  u  — ►  go.  (4.11) 

In  order  to  investigate  the  limiting  process  |v|  _»  oo  we  restrict 
M  by 


M  |  >  1. 


(4.12) 


Thanks  to  this  inequality  we  are  assured  that  for  sufficiently  large 
values  of  J  v  J  the  inequalities  (4. Id)  are  satisfied.  We  obtain 

.2 


Fq  — *  iw  a<)  po 


l(M  -1)1  MZ-1  J  l(M  -1)5  J 


2  2 

\  -  -  ;  F  —  0.  m  «  2,3;  F*  -  0,  m  «  1,2,3  (4.13a) 

1  c  ..c  ,  m  m 

M  -1 


for  v  I  —  oo. 


These  relations  lead  to 


Ap-{4%a,p“ 

—  w^*  0  for  j  v  |  -►  oo  • 


2  2 
ao  Po  M 

IT 


(4.13b) 


From  the  relation  (4.13b)  we  obtain  immediately 


^  -  (4*  -  ^  »)  -  0  for  |  M  |  — *■  co ,  |v|-*oo,  (4.13c) 
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Ap  -  Ap  —  Ofor|M|  —  oo,  j v  J  —  oo,  R  —  oo.  (4.13d) 

In  case  the  inequality  (4.12)  does  not  hold,  the  restrictions  (4.  Id) 

would  be  violated  sooner  or  later  as  jvj  tends  to  infinity,  no  matter 

how  large  the  value  of  u  is.  Then  the  expansion  (4.1a)  is  no  longer 

valid.  This  case  will  be  dealt  with  shortly  in  section  5. 

2 

It  should  be  emphasized  that  the  coefficient  -  of  the  cor¬ 

rection  term  -  w,  which  appears  especially  in  eqs.  (4.7b),  (4.10b), 
and  (4.13c)  is  independent  of  the  parameters  M,  «  ,  v,  and  n. 

It  is  almost  superfluous  to  quote  the  following  remark:  The 

asymptotic  expansions  in  section  3  were  carried  out  up  to  the  term 

.3 

of  the  order  £  .  In  accordance  with  this  the  expansions  (4.1a)  and 

(4.6 )  are  extended  up  to  the  term  of  the  order  R"  .  In  case  the  expan¬ 
sions  would  be  carried  out  to  higher  order  terms,  one  would  obtain 
expansions  for  the  aerodynamic  pressure  analogous  to  (4.1a)  and  (4.6). 
Formulas  (4.  3c),  (4.3d),  (4.4e),  (4.4f),  (4.  7)  to  (4.13)  would  be  valid 
too,  after  the  range  of  the  index  m  is  extended  appropriately. 

From  eqs.  (4.1a),  (3.4),  and  (3.5)  one  learns  that  the  coeffi- 
cients  of  (l/R)  ,  (l/R)  in  the  asymptotic  expansion  for  Ap  depend 
upon  n^.  (The  same  is  true  for  the  coefficients  of  higher  powers  of 
l/R,  if  the  asymptotic  expansion  would  be  carried  out  to  higher  terms, 
see  section  3.)  For  the  investigation  of  the  flutter  of  cylindrical 
shells  it  is  often  necessary  to  consider  values  of  n  up  to  the  order 
of  20,  see  refs.  1  and  6.  For  these  large  values  of  n  the  influence 
of  the  terms  of  the  order  (l/R)  and  (l/R)  in  (4.1a)  and  (4.6 )  can 
become  quite  significant. 
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5.  Examination  of  the  Remaining  Cases. 

In  section  4  only  the  case  |  j  >  1,  j  j  >  1  has  been 
considered,  where  the  linear  piston  theory  expression  proved  to  be 
a  first  order  approximation  for  the  aerodynamic  pressure  under  in¬ 
vestigation.  In  the  remaining  case,  i.e. 


a.)  |Mj  <1;  |M2|<1; 
c.)  |  MjJ  >  1;  jM2|<l 


b.)  |mJ<1;  j  M2  J  >  1; 

(5.1) 


the  formulas  (2.11),  (2.12)  and  similar  ones  have  to  be  applied,  and 
besides  the  expansions  (3.4)  and  (3.5)  also  the  asymptotic  expansion 
(3.6)  has  to  be  used.  In  none  of  these  remaining  cases  (5.1)  the 
linear  piston  theory  expression  (4.2)  can  be  considered  as  a  first 
order  approximation  of  the  aerodynamic  pressure  under  investigation. 
Let  w  ,  the  cylinder  deformation,  again  be  given  by  eq.  (2.  4).  Then 
one  obtains  for  instance: 


Ap  -  Ap  s 

«  ao  Po  ji  (  iM  +  2M  ^  ^  +  . . .  )  -  (5.2a) 

-v(-^*  “2+  + m?1  *  •••)  W+°<R'1>} 

for 

|Mj  |  c  !;  |M2  j  <  1;  v>  0  (5.2b) 


and 
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** 

Ap - Ap  m 


i  a02Po{-v(M12+iM14+...)w^  iv(M-  ^-+ i  .  )  w  - 

^  (5  •  3a) 

-i(Mf+iM14+...)  (  •••)  |j  + 

2 

+0(R-1)j- 


for 


|  Mj  1  <  1;  M2  >  Ij  v  >  0 


(5.3b) 
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6.  Conclusions. 

In  the  preceding  sections  the  aerodynamic  pressure  acting  on 
a  circular  cylindrical  shell  of  infinite  length  in  an  airflow  parallel  to 
the  cylinder  axis  is  investigated.  The  shell  is  deformed  by  a  har¬ 
monically  oscillating  standing  wave  of  a  sinusoidal  pattern.  Based 
upon  the  exact  solutions,  furnished  by  Leonhard  and  Hedgepeth,  asymp¬ 
totic  expansions  are  developed  for  the  aerodynamic  pressure.  These 
expansions  allow  an  investigation  of  the  accuracy  of  the  linear  piston 
theory  approximation  (4.2)  when  applied  to  cylindrical  shells.  Fur¬ 
thermore  improved  approximations  can  be  obtained  from  these 
asymptotic  expansions. 

In  case  the  absolute  values  of  and  (see  eqs.  (2.2) 

and  (2.  5) )  are  larger  than  one,  |  MjJ  >  1;  |  |  >  1,  then  the  linear 

piston  theory  expression  Ap  can  be  considered  as  a  first  order  ap¬ 
proximation  for  the  aerodynamic  pressure  under  investigation.  Never- 
theless  the  replacement  of  Ap  by  the  approximation 

***  **  aoPo  f  x,  dw ,  9w  ao  1  .t, 

Ap  *  Ap  -  w  =  a0p0  |  ao  M^+ ^  -  ^-wj-  (6.1) 

is  suggested  as  a  first  step  of  an  improvement  for  the  application  to 

2 

cylindrical  bodies.  The  coefficient  -  ^p0  of  the  additional  term  in 
(6.1)  is  independent  of  the  parameters  M,v,w  and  n.  Due  to  the 
independence  of  v  the  approximation  formula  (6.1)  can  be  applied  to 
more  arbitrary  oscillatory  shell  deformations  of  the  form 

w(x,8,t)  «  cos(n6)  w(x)  eiwt  ,  (6*2) 

where  w(x)  and  its  first  derivative  are  representable  as  Fourier 


-25- 


aeries  (provided  that  for  all  individual  terms  the  above  mentioned 
inequalities  for  j  M^  |  and  J  M2  J  are  valid)  without  a  constant  term. 
The  last  remark  reflects  the  fact  that  v  *  0  was  excluded  from  our 
investigations. 

For  large  values  of  n  especially  the  correction  terms  of 
the  order  (l/R)^  and  (l/R)^  in  eqs.  (4. la)  and  (4. 6)  can  become 
quite  significant.  For  flutter  investigations  it  is  often  necessary  to 
consider  values  of  n  up  to  the  order  of  20. 

In  the  following  cases  a. )  J  M1  ]  <  1;  j  M2  j  <  1;  b. )  |  Mj  |  <  lj 
I  ^*2  1  >  c* )  |  Mj  J  >  1;  j  M2  j  <  1  (see  (5.1))  the  linear  piston  theory 
expression  (4.2)  can  no  longer  be  considered  as  a  first  order  approxi¬ 
mation  for  the  aerodynamic  pressure,  as  is  demonstrated  in  section  5. 
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